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Abstract. We analyze the effects of different forms of the sound-velocity 
function c s (T) on the hydrodynamic evolution of matter formed in the central 
region of relativistic heavy-ion collisions. At high temperatures (above the 
critical temperature T c ) the sound velocity is calculated from the recent lattice 
simulations of QCD, while in the low temperature region it is obtained from 
the hadron gas model. In the intermediate region we use different interpolations 
characterized by the values of the sound velocity at the local maximum (at T = 
0.4 T c ) and local minimum (at T = T c ). In all considered cases the temperature 
dependent sound velocity functions yield the entropy density, which is consistent 
with the lattice QCD simulations at high temperature. Our calculations show 
that the presence of a distinct minimum of the sound velocity leads to a very 
long (~ 20 fm/c) evolution time of the system, which is not compatible with 
the recent estimates based on the HBT interferometry. Hence, we conclude that 
the hydrodynamic description is favored in the case where the cross-over phase 
transition renders the smooth sound velocity function with a possible shallow 
minimum at T c . 
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The recently analysed RHIC data suggests that matter formed in relativistic 
heavy-ion collisions behaves like a perfect fluid |T|. This hypothesis triggers new 
developments in relativistic hydrodynamics I21I3HlISlISlC7lIHlia HDl ttH H21 UHl H3J - 
In this paper we continue the studies of Ref. [TS] and discuss the effects of different 
forms of the sound- velocity function c s (T) on the hydrodynamic evolution of matter 
formed in the central region of relativistic heavy-ion collisions. 
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Figure 1. The sound velocity squared, c 2 s (T), shown as a function of the 
temperature at zero baryon density. The left panel shows the results of the lattice 
simulations of QCD (solid line) and the ideal hadron gas model [17] (dashed 
line). The simplest interpolation between the two calculations is marked by a 
thick solid line. The right panel presents the three cases of the function c s (T) 
considered in our calculation. The case I agrees with the interpolation shown in 
the left panel. A deeper local minimum at T c is introduced in the case II (dashed 
line) and III (long-dashed line) to model different forms of the phase transition. 
The sound velocity at the critical temperature is 25% (case II) and 50% (case III) 
lower in comparison to the case I. A corresponding increase of the sound velocity 
around T sa 0.4 T c in the cases II and III is required to have the same entropy 
density at high temperatures. The increase of the sound velocity in this region 
may account for pourly known repulsive van der Waals forces in the hadron gas. 

The left panel of Fig. [T] shows the sound velocity function c 2 s (T) obtained from the 
lattice simulations of QCD [16] (solid line) and the ideal hadron gas model [17] (dashed 
line) . The lattice results were obtained for physical masses of the light quarks and the 
strange quark. One can observe that the sound velocity is below the limiting value 
c 2 s = | expected for a non-interacting gas of massless particles. On the other hand, in 
the low temperature part the hadron-gas calculation is used which includes full mass 
spectrum of resonances [18j . It is interesting to check that for very small temperatures 
(T << to?,-) the sound velocity is determined by the properties of a massive pion gas 
giving c s {T) = y/T/m^, where is the mass of the pion. The lattice and hadron- 
gas calculations overlap in the region close to the critical temperature T c = 170 MeV 
and can be naturally joined by an interpolating function (thick solid line). A smooth 
function like that would result in a cross-over phase transition. A different way to 
connect the two results is to introduce a local minimum at the critical temperature. 
As the value of the c s decreases, the phase transition more and more resembles the 
first order phase transition where, by definition, the sound velocity drops to zero at 
T = T C . 

Since the lattice QCD simulations are not very much reliable in the region 
T < T c and the hadron-gas calculations are meaningless in the region T > T c we 
decided to consider different interpolations between the lattice and hadron-gas results 




Figure 2. The temperature dependence of the entropy density and energy 
density, the two upper panels, as well as the energy density dependence of the 
pressure and sound velocity, the two lower panels. One can observe that the 
deeper is the minimum of the sound velocity function, the steeper is the increase 
of the entropy density and the energy density. 



in the vicinity of the phase transition. In this respect our study may be regarded 
as supplementary to Ref. [T^] where the impact of different equations of state on 
the particle spectra and V2 was studied. The three considered by us interpolations 
are denoted below as the case I, II and III, for details see Fig. [TJ All considered 
functions are consistent with the lattice simulations of QCD for the temperatures 
T > 1.15 T c and with the ideal hadron gas for the T < 0.15 T c . In the temperature 
range 0.85 T c < T < 1.15 T c the three cases differ by the depth of the local minimum. 
The value of the sound velocity at the critical temperature T c equals 0.37 in the case 
I. In the case II it is decreased by 25% compared to the case I (c s (T c ) = 0.28) and in 
the case III it is decreased by 50% (c s (T c ) — 0.19). In order to be able to compare 
the solutions obtained for such three different cases, it is crucial to have the same 
entropy density at high temperatures, where the lattice data are reliable. To achieve 
this requirement we use the following thermodynamic relation connecting the entropy 
density s(T) with the sound velocity for zero baryon chemical potential 

T 

dT 



s(T) = s(T ) exp 



To 



T'ci(T') 



(1) 



One may conclude from Eq. |T]) that the decrease of the sound velocity results 
in the relative increase of the entropy density in the high temperature region. To 
avoid the entropy increase we increase the values of the sound velocity in the region 
0.15 T c < T < 0.85 T c . This kind of modification may be interpreted as an effective 
parametrization of the repulsive van der Waals forces in the hadron gas. In Fig. [2] we 
present several thermodynamic functions calculated for the 3 considered sound velocity 
functions. The upper left plot shows the functions s(T)/T 3 where, as required, all the 
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lines converge at high temperatures. All presented sound velocity functions have been 
also checked to fulfill the condition against the formation of shock waves [30] HI] 

A detailed explanation of the method to solve the relativistic hydrodynamic 
equations with the temperature dependent sound velocity may be found in Ref. [20] 
and in our later developments of this approach in [17l [22], see also [23]. When 
considering RHIC experiments in the midrapidity region one can assume with good 
justification boost-invariance and zero baryon chemical potential [24] [25] . With this 
simplification the only independent thermodynamic parameter is the temperature and 
the equations of relativistic hydrodynamics have the following form 

u^(Tu v )=d v T, d„(sO = 0, (2) 

where T is the temperature, s(T) is the entropy density, and u M = 7(l,v) is the 
hydrodynamic four-velocity. Further transformations include introduction of the 
temperature dependent sound velocity c s (T) and the potential $ (T) according to 
the formulas 

2 DP s dT ^ d\nT 

C ^7fe=T97' <» = — = (3) 

The relativistic boost-invariant hydrodynamic equations with cylindrical symmetry 
may be expressed in the compact form with the help of the auxiliary functions A + 
and A_ defined as 

A ± = <f>±a, (4) 

where a is the transverse rapidity of the fluid element, v r — tanha. The longitudinal 
velocity in the considered case has the well known form v z — z/t jJS] and Eqs. (|2J) 
read 

dA± (t, r) + v r ±c s dA± (t,r) + c s fv^ + 1\ = Q 



dt 1 ± v r c s dr 1 ± v r c s \ r t r 

All interesting physical observables can be calculated from the auxiliary functions A± 
with the explicit formulas given in Refs. [i"5 ] IT7 l [22], 

For symmetry reasons, both the velocity field v r and the temperature gradient 
dT/dr should vanish at r = 0. This condition is achieved if the functions A + and A_ 
are initially determined by a single function A(r) according to the prescription [20] 

A+(t = t ,r)=A(r), A_(t = t„,r) = A{-r). (6) 

Our main physical assumption about the initial state is that the initial temperature 
profile is connected with the nucleon- nucleus thickness function Ta(t) by the following 
equation 

Ta(t) 
'Ta(O) 

Here T s (s) is the inverse function to the entropy density function s(T) and the 
parameter so is the initial entropy at the center of the system. The idea to use Eq. 
([7]) follows from the assumption that the initially produced entropy density s(to,r) 
is proportional to the density of wounded nucleons at a distance r from the collision 
center. We use the value s = 70.5 fm -3 which yields T(t = t ,r = 0) = 2T C . We 
note that the functions s(T) and T s (s) (evaluated for the cases I, II, and III) agree 
if the temperature or entropy is sufficiently large. Hence, in the three considered 
cases the initial temperature profiles are practically the same in the center of the 
system. The small differences appear however if we consider larger values of r where 



T{t ,r) =T s [s(t ,r)} =T S 



So; 



(7) 
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Figure 3. Isotherms describing the hydrodynamic evolution for the case I; sound 
velocity with a shallow minimum where c s (T c ) = 0.37; case II sound velocity 
with a moderate minimum, c s (T c ) = 0.28 and case III with a deep minimum, 
c s (T c ) = 0.19. The numbers denote the values of the temperature in the units of 
the critical temperature T c = 170 MeV. 



the temperature and entropy drops down. We recall that the thickness function is 
defined by the equation 

T A (r) = 2 J dzp(Vr*+z 2 ), (8) 

where the function p(r) is the nuclear density profile given by the Woods- 
Saxon function with a conventional choice of the parameters: po — 0.17fm -3 , 
r Q = (1.12A 1 / 3 - O.86A- 1 / 3 ) fm, a = 0.54 fin, A = 197. We note that the initial 
condition ([Jj may be included in the initial form of the function A(r) 

(9) 

In Fig. [3] we show the isotherms describing the hydrodynamic evolution of the 
system in the cases I, II and III. The numbers at the isotherms give the values of the 
temperature in the units of the critical temperature T c . In the case I, the upper panel, 
we observe that the center of the system cools down to 0.8 T c after the evolution 
time of about 15 fm. We note that the hydrodynamic description should be replaced 
(around T ~ 0.8 T c ) by the model describing hadronic rescattering (see, e.g., Ref. [IT] ) 
whose presence additionally increases the lifetime of the system. Another option is to 
assume that freeze-out happens at high temperature (see Refs. [371 HE]) where many 



A(t = t ,r) = <S>t<T : 



so 



Ta(t) 
Ta(0)_ 
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physical observables were successfully reproduced under the assumption of a universal 
freeze-out taking place at the temperature of 165 MeV). In the latter case the lifetime 
of the system may be identified simply with the time when the system passes the 
phase transition. In the discussed case this time is of about 10 fm. 

In the lower left panel of Fig. [3] we show the isotherms describing the 
hydrodynamic evolution of the system with the sound velocity II. The initial entropy 
is exactly the same as in the case I. One can notice that the central temperature does 
not drop below 0.8 T c before 20 fm. Clearly, the dip in the sound velocity causes a 
dramatic increase of the lifetime of the system. The most striking situation is presented 
in the lower right part of Fig. [3]where the sound velocity with the deepest minimum is 
considered. One can notice that the system does not pass the phase transition before 
the considered evolution time of 20 fm. Contrary, even at that time the system has 
the tendency to expand more, the effect indicated by the shapes of the isotherms. 

The observation of the effects related with the possible presence of the softest 
point of the QCD equation of state triggered the ideas about the long-living states 
formed in ultra-relativistic heavy-ion collisions [29l [30j [31] . In view of the RHIC data 
indicating a short lifetime of the system we may analyze the problem of the lifetime 
in the reverse order, asking the question how soft the equation of the state may be 
to allow for a hydrodynamic description consistent with the HBT results. Our results 
show that the smooth behavior of the sound velocity without any distinct minimum 
is favored if we demand the short hydrodynamic evolution time of about 10 fm. 



U. W. Heinz, |nucl-th/051205l] 

D. Teaney, Phys. Rev., C68 (2003) 034913. 

T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey, Y. Nara, Phys. Lett., B636 (2006) 299. 
Y. Hama, et al, Nucl. Phys., A774 (2006) 169. 

K. J. Eskola, H. Honkanen, H. Niemi, P. V. Ruuskanen, S. S. Rasanen, Phys. Rev., C72 (2005) 
044904. 

U. W. Heinz, H. Song, A. K. C haudhuri, Phys. R ev., C73 (2006) 034904. 
P. Huovinen, P. V. Ruuskanen, |nucl-th /0605008 

R. Andrade, F. Grassi, Y. Hama, T. Kodama, J. Socolowski, O., Phys. Rev. Lett., 97 (2006) 
202302. 

R. Baier, P. Romatschke, U. A. Wiedemann, P hys. Rev., C73 ( 2006) 064903. 
T. Koide, G. S. Denicol, P. Mota, T. Kodama, |hep-ph/0609117| 
C. Nonaka, S. A. Bass, Phys. Rev., C75 (2007) 014902. 

L. M. Satarov, I. N. Mishustin, A. V. Merdeev, H. Stoecker, hep-ph/0611099 
T. Hirano, U. W. Heinz, D. Kharzeev, R. La cey, Y. Nara, [n ucl-th/0701075 
A. Bialas, M. Chojnacki, and W. Florkowski, larXiv:0708. 10761 |nucl-th|. 
M. Chojnacki, W. Florkowski, |Wl-th/0702030") 

Y. Aoki, Z. Fodor, S. D. Katz, K. K. Szabo, JHEP, 01 (2006) 089, |hep-lat705 10084 
M. Chojnacki, W. Florkowski, T. Csorgo, Phys. Rev., C71 (2005) 044902. 
W. Broniowski, W. Florkowski, Phys. Lett., B490 (2000) 223. 
P. Huovinen, Nucl. Phys. A761 (2005) 296. 

G. Baym, B. L. Friman, J. P. Blaizot, M. Soyeur, W. Czyz, Nucl. Phys., A407 (1983) 541. 

J. P. Blaizot, J.-Y. Ollitrault, Phys. Lett., B191 (1987) 21. 

M. Chojnacki, W. Florkowski, Phys. Rev., C74 (2006) 034905. 

A. Dyrek and W. Florkowski, Acta Phys. Pol. B15 (1984) 653. 

W. Florkowski, W. Broniowski, M. Michalec, Acta Phys. Polon., B33 (2002) 761. 

G. Torrieri, et al, Comput. Phys. Commun., 167 (2005) 229. 

J. D. Bjorken, Phys. Rev., D27 (1983) 140. 

W. Broniowski, W. Florkowski, Phys. Rev. Lett., 87 (2001) 272302. 

W. Broniowski, A. Baran, W. Florkowski, Acta Phys. Polon., B33 (2002) 4235. 

C. M. Hung, E. V. Shuryak, Phys. Rev. Lett., 75 (1995) 4003. 

D. H. Rischke, M. Gyulassy, Nucl. Phys., A597 (1996) 701; 
U. W. Heinz, J. Phys. Conf. Ser., 50 (2006) 230. 



